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The boson-fermion (BF) model has been proposed to describe super- 
conductivity in short coherence length superconductors. In this work we 
use it to study impurity induced inhomogeneities in high temperature su- 
perconductors (HTS) as found in numerous scanning tunnelling measure- 
ments. The model has been formulated in real space and solved with help of 
Bogolubov-de Gennes approach. Disorder in the boson or fermion subsys- 
tem is directly coupled to the superconducting order parameter and leads 
to severe changes of superconducting properties like local order parameter 
and density of states. We present the results for many impurities randomly 
distributed over otherwise clean and periodic two-dimensional square lat- 
tice. 

PACS numbers: 71.10.Fd; 74.20.-z; 74.81.-g 



1. Introduction 

The boson-fermion model assumes the existence in the system of two 
kinds of interacting quantum objects PQ consisting of (usually) local and 
immobile bosonic particles (doubly charged preformed pairs) interacting 
with mobile fermions (electrons). It is the boson-fermion scattering which 
induces the condensation transition in bosonic and superconducting tran- 
sition in fermionic subsystem |2j. In real materials, as e.g. doped HTS, 
one expects both bosonic and fermionic parameters to depend on the local 
environment and thus to be random variables. 

Here we present results of calculations of such important parameters 
as local value of the gap parameter Aj and local density of states Ni(E) 
at lattice site %. A, and Ni(E) are subject to measurements with help of 
scanning tunnelling microscope (STM). The STM spectra of many high 
temperature superconductors [3], e.g. Bi2Sr2CaCu20s+,5, indeed show the 
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intrinsic inhomogeneities of A, and Ni(E) extending over the scales of few 
nanometres. 

Even though the microscopic origin of these inhomogeneities remains 
unclear some experiments show positive correlation between location of the 
oxygen dopant and the value of local order parameters. Based on this sce- 
nario theoretical approaches have been proposed in which authors calculate 
local properties of the system by means of real space Bogolubov-de Gennes 
approach, using Hubbard |4j, t-J [5] or similar models with disorder in nor- 
mal (density) or anomalous (Cooper) channel [6]. To account for possible 
sources of inhomogeneities we use the boson-fermion model [I] in the form 

^ BF = E >u<i<> + E (5 mp - m) + (i) 

+ + 8Ef - 2/x) b% + Y J 9 f ('>,'■ < r + &44) , 

i hj 

where i and j denote lattice sites of the square lattice, c\ a (ci, CT ) stand for 
creation (annihilation) operator of fermion at the site i with spin a and b\ 
and bi are creation and annihilation operators of hard-core bosons at the 
site i. fj, denotes chemical potential of the system and Uj are hopping in- 
tegrals. For numerical calculations we assume hopping integrals to nearest 
and next nearest neighbours ti,t% to be different from zero. As noted earlier 
the presence of oxygen impurities in Bi family of HTS on one hand is the 
source of free carriers (holes) and at the same time introduces disorder into 
the materials. Undoped system is insulating. We assume that this disorder 
can be modelled by random on-site energies V^ imp in fermion subsystem and 
random hard-core boson energies bEf [3 [8]. To account for the short coher- 
ence length and d-wave symmetry of the superconducting order parameter 
we assume that the boson-fermion coupling takes on non-zero value for 
nearest neighbour sites < i, j > only and is equal +g if j = i ± x and -g if 
j = i± y. 

Application of standard Hartree-Fock-Bogolubov decoupling leads to 
H BF = H B + H F , where H B is single site bosonic hamiltonian with param- 
eters depending on the fermionic order parameter. H F is the Hamiltonian 
describing fermionic subsystem [8] . Standard statistical approach allows the 
exact solution of the boson Hamiltonian H B . One finds 

<5,> = -^tanh(^). (3) 
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Here H = \j ( ~ +<5 f' : — ) + Ixd 2 , Xi = < 6 U%,T > and A *i = 

<i> 

< bj >. It is Ajj which couples two subsystems as is evident from 
equation fl3J) below. 

The fermion part has the standard BCS-structure and we diagonalise it 
by the Bogolubov-Valatin transformation, which yields the following Bogolubov- 
de Gennes equations [8 J 

( u \ 




We solve numerically the above Bogolubov-deGennes equations in real 
space iteratively to get self-consistent values for all Aj,- with high accuracy. 
Even though our model and calculations are motivated by the experimental 
results obtained mainly for Bi family high temperature superconductors we 
do not try to model the actual experimental situation here. Instead we shall 
present the results on local properties of disordered boson-fermion model. 
In this respect the paper is an extension of the previous CPA calculations 
[8J of the same model. 



2. Results and discussion 

Before we start the discussion of the results let us remind that it is the 
position of the boson level with respect to the Fermi level Eb which tells us 
if the system will be superconducting at all. Due to phase space restrictions 
at low temperatures the scattering of preformed bosons is allowed only, if 
the bosonic level is close enough to the Fermi level. 

To account for disorder in the system the boson energies Ef = E B +5Ef 
were randomly and uniformly distributed around level E B over the interval 
AE B . Physically this corresponds to fluctuating pairing potential [6], which 
in BF model arises in a very natural way. 

Due to the above mentioned property of the model the magnitude of the 
local order parameter < bi > (not shown) obtained as a solution of equation 
@ is correlated with the position of the impurity provided we identify the 
dopant oxygens as e.g. those in Bi2Sr2CaCu20g + 5 with the pairing centres. 
From the same calculations it can also be seen that both order parameters 
< bi > and Xi take on large values in the same regions of space. The direct 
proportionality of both order parameters is a general property of the model 
[2], which can also be inferred from the mean- field equation ([3]). 

In figure (1) we show the effect of disorder in bosonic (random Ef levels 
- left panel) and fermionic (random V^ mp - right panel) subsystems on the 
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Fig. 1. Average values of | < 6, > | and as a function of disorder in the 
boson (left panel) and fcrmion (right panel) subsystem. Inset to left panel: The 
temperature dependence of of average < hi > | compared with similar dependence 
at the sites at which the value of the gap is larger /smaller than average. The results 
were obtained for t x =t = 1.0, t 2 = -0.3i, [i = O.Ot, g = 0.8t E B = -0.2t. 



average values of the order parameters | < bi > | and \xi\- In this figure the 
clean system (characterised by Eb = — 0.2£) would correspond to the ex- 
perimental situation with all sites contributing to boson-fermion scattering. 
With increase of AE B the number of sites at which this scattering takes 
place diminishes and the order parameter is strongly reduced. 

It is interesting to note that despite the d-wave character of supercon- 
ductivity the disorder acting on the fermionic subsystem relatively weakly 
affects both order parameters (figure (1), right panel). More severe changes 
of both | < bi > | and \xi\ are observed for very large disorder. The results 
presented in inset to figure (1) show that the whole system is characterised 
by the unique transition temperature. This result is correct from thermo- 
dynamical point of view. At the same time it indicates that for the system 
studied here there exist a coupling between regions of large and small pair- 
ing potential, preventing the appearance of isolated regions with different 
values of the transition temperatures, even though the gaps near T c do vary 
a lot as is also visible from figure (2) - right panel. Indeed rough estimation 
of the coherence length leads to values of order 1/3 to 1/2 of the system 
size. This shows that all domains are coupled. In figure (2) we plot the 
maps of bosonic order parameter | < 6j > | of our 23x23 large system 
at different temperatures. They are shown in gray scale for two values of 
temperature - T « T c (left panel) and for T s=s 0.9T C (right panel). At 
temperatures close to T c the regions of large gap seem to be surrounded by 
the non-superconducting material. Close inspection shows that the whole 
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Fig. 2. Maps of amplitude of bosonic order parameter | < 6» > | for different values 
of temperature: T = O.Oi (left panel), T = O.li w 0.9T C (right panel). Other 
parameters the same as in figure (1) 



system is superconducting but most of it is characterised by very low value 
of the gap, while large gap persists in few regions of the sample. We do 
not address the pseudogap problem here, as our mean-field approach is not 
suitable in that region of parameter space. 

In conclusion, we studied the properties of disordered boson-fermion 
model allowing for two different disorder sources. The randomness in boson 
energies strongly degrades superconductivity, while disorder acting on elec- 
tron subsystem y. imp has much weaker effect on superconducting state. The 
local values of fermionic order parameter positively correlate with positions 
of "pairing impurities" characterised by small values of Eg. Even strongly 
disordered system possesses single superconducting transition temperature 
despite large fluctuations of the gap magnitude. 
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